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ABSTRACT 


Two  current  vector  time  series  obtained  in  the  Mississippi  Bight  exhibit  clock- 
wise polarized  currents  of  near-inertial  frequency  that  are  closely  associated  with 
shifting  winds.  Because  of  the  closeness  of  the  local  inertial  period  and  the 
diurnal  tidal  period,  it  is  difficult  at  first  glance  to  determine  the  true  nature 
of  the  observed  rotary  currents.  However,  complex  demodulation  at  the  Inertial 
frequency  reveals  a strong  signal  accompanying  wind  shifts  that  are  usually  asso- 
ciated with  the  passage  of  atmospheric  fronts. 

Spectral  analysis  for  clockwise  and  counterclockwise  frequencies  indicates  a 
highly  energetic  peak  in  the  inertial-diurnal  frequency  band  for  the  clockwise 
spectrum.  The  rotary  coefficient  computed  from  the  autospectra  and  quadrature  spec- 
trum of  the  vector  components  gives  C[^  > 0.  9 in  the  vicinity  of  the  inertial-diurnal 
frequency  band. 

A model  using  wind  stress  as  a forcing  function  is  highly  effective  in  repro- 
ducing sinusoidal  oscillations  seen  in  the  observed  current.  These  oscillations 
occur  in  conjunction  with  shifts  in  the  wind  direction. 

Because  of  the  close  association  of  the  near- inertial  oscillations  with  local 
wind  effects,  it  is  concluded  that  inertial  currents  are  locally  induced  by  wind 
stress.  Furthermore,  wind  stress  not  only  initiates  the  rotary  currents  but  is 
also  highly  effective  in  destroying  n. 
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I.  INTRODUCTION 

Inertial-period  oscillation  of  horizontal  currents  is  a phenomenon  that  has 
been  observed  at  a variety  of  depths  and  locations  on  the  globe.  The  inertial 
period,  Tj , is  equal  to  2Il/f,  where  f is  the  Coriolis  parameter,  defined  as  f = 
2f2sin(().  The  quantity  SI  is  the  angular  velocity  of  the  earth's  rotation  and  is 
the  geographical  latitude.  The  inertial  period  varies  from  a value  of  12  hours  at 
the  poles,  to  24  hours  at  30°  latitude,  to  infinity  at  the  equator. 

Under  ideal  conditions,  inertial  oscillations  in  a horizontal  current  field 
occur  under  a free-flow  regime  when  no  forces  act  on  the  fluid  other  than  the 
Coriolis.  The  effect  of  Coriolis  force  is  to  deflect  tlie  motion  to  the  right  in 
the  northern  hemisphere,  and  therefore  its  only  influence  is  in  changing  the  flow 
direction  and  not  the  magnitude.  Under  these  conditions,  water  flows  in  horizontal 
circular  motion,  clockwise  in  the  northern  hemisphere  and  counterclockwise  in  the 
southern  hemisphere  with  period  of  rotation  Tq. 

The  horizontal  equations  of  motion  (Neumann,  1968;  Neumann  and  Pierson,  1966) 
for  the  case  of  accelerated  currents  are 


du  _ 
dt 


fv 


3P 


dv 

dt 


= -fu  - F 


y 


(1) 


where  u is  the  speed  in  the  x direction,  v is  the  speed  in  the  y direction, 
d/dt  = 3/3t  + u 3/3x  + v 3/3y  is  the  so-called  "material  derivative,"  f is  the 
Coriolis  parameter,  Fjj  and  Fy  are  the  frictional  forces  along  the  x and  y direc- 
tions, respectively,  a is  the  specific  volume,  and  P is  the  pressure.  If  the 
Guldberg-Mohn  assumption  is  used  to  describe  the  frictional  forces,  then  F^^  = ku 
and  Fy  - kv,  where  k is  a proportionality  constant. 

Multiplying  the  first  equation  by  u and  the  second  by  v and  adding  yields 


du 

" ^ ^ ' dt 


^ _ 1 dc  _ / 3P  , ,2 

■ 2 di “ air 


(2) 
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where  c = u + v*^. 

For  a homogeneous  ocean,  the  pressure  gradients  in  equation  2 can  result  only 
from  a slope  of  the  sea  surface.  Although  the  slope  of  the  sea  surface  can  be  a 
significant  driving  force  in  coastal  regions,  the  site  of  this  study  seems  to  be 
sufficiently  far  from  the  nearest  coast  (about  25  km)  that  slope  effects  are  negli- 
gible, at  least  in  the  diurnal  frequency  band.  The  model  used  in  this  report 
assumes  that  no  pressure  effects  were  present,  yet  it  gives  excellent  correspon- 
dence between  observed  and  modeled  currents  (see  Section  II-E). 
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On  the  other  hand,  sea  surface  slopes  represent  potential  energy  tliat  may  later 
be  converted  to  kinetic  energy  of  the  inertial  currents.  For  example,  if  a wind 
stress  that  sets  up  and  maintains  a slope  is  suddenly  removed,  the  slope  may  rapidly 
disappear  and  the  stored  energy  will  be  converted  to  current  motion.  Such  current 
motion  is  acted  upon  mainly  by  friction  and  Coriolis  force.  The  result  is  a damped 
current  rotating  at  inertial  frequency.  Therefore,  if  we  require  the  horizontal 
pressure  gradient  to  go  to  zero,  equation  2 becomes 


c dt 


-k  . 


(3) 


Integration  of  this  equation  yields 


c = c exp  (-kt)  (A) 

o 

where  Cq  is  the  current  speed  at  t = 0.  For  k = 0,  equation  3 is  the  equation  for 
a pure  frictionless  inertial  current.  For  k 0,  equation  3 represents  a damped 
inertial  current.  Equation  4 indicates  that  decay  of  current  magnitude  is  an 
exponential  function  of  k.  The  value  of  k can  vary  from  place  to  place  in  the 
ocean;  typical  values  of  10”^  to  10“^  sec"l  have  been  reported  by  Defant;  0.7  x 
10“5  sec”3,  by  Gustafson  and  Kullenberg  for  the  Baltic  Sea;  and  3 x 10"^  sec”^,  by 
Neumann  for  the  upper  thousand  meters  of  the  ocean  (Neumann,  1968). 

If  we  multiply  the  first  part  of  equation  1 by  v and  the  second  by  u and  then 
subtract,  we  obtain,  for  no  horizontal  pressure  gradient. 


= f(u2  + v2) 

dt  dt 


Now,  d/dt  (u/v)  = (vdu/dt  - udv/dt)/v2. 


Therefore, 


(5) 


2 f 

dt  V 2 
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In  Cartesian  coordinates  u/v  = cotB  and  v“^  = c sin  B,  where  B is  the  angle  between 
the  x-axis  and  the  current  direction.  Thus, 


f 

2 

sin  B 


or 


= -f  . (7) 

dt 

Equation  7 states  that  the  time  rate  of  change  of  the  direction  of  an  inertial  cur- 
rent is  a constant  that  is  the  negative  Coriolis  parameter.  In  the  northern 
hemisphere  iji  > 0,  f > 0,  and  therefore  dB/dt  < 0,  l.e.,  clockwise  motion. 

The  present  investigation  makes  use  of  a 32-day  (28  January  to  28  February 
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1974)  moored  current  meter  record  obtained  in  the  Mississippi  Bight  by  members  of 
the  Coastal  Studies  Institute  as  part  of  the  I.ouisiana  Offshore  Oil  I'ort  (LOOP; 
project  (Wiseman  et  al.,  197b).  A second  record,  comprising  39  days,  is  also  dis- 
cussed, but  analysis  was  limited  because  of  lack  of  meteorological  data. 

The  current  meter  was  moored  south  of  Grand  Isle,  Louisiana,  at  latitude 
28°55'N,  longitude  89°47'W, and  6.4  meters  below  the  surf  ace  in  water  33  meters 
deep.  Wind  data  was  obtained  from  an  anemometer  placed  on  a Chevron  rig  located 
about  15  km  northwest  of  the  site  of  the  current  meter  at  latitude  28‘’01'N, 
longitude  90°06'W.  (See  Figure  1 for  locations.) 

It  is  the  purpose  of  this  study  to  show  that  inertial-period  motions  in  the 
study  area  are  excited  by  strong  winds,  in  particular  those  winds  that  accompany 
atmospheric  frontal  passage.  The  major  difficulty  in  quantitatively  discerning 
I these  motions  from  others  arises  from  the  fact  that  the  inertial  period  at  the 

I latitude  of  the  study  site  is  24.76  hours,  which  corresponds  approximately  to  the 

• diurnal  tidal  period.  This  problem  has  plagued  a number  of  researchers  who  have 

j worked  at  latitudes  near  30°.  Webster  (1968),  in  a comprehensive  review,  listed 

24  significant  studies  that  involve  inertial-period  oscillations  spanning  the  years 
I 1931  to  1966.  He  noted  that  in  the  older  literature  it  is  often  difficult  to 

I determine  whether  what  the  researchers  reported  were  indeed  inertia]  motions.  On 

1 the  other  hand,  certain  investigators  may  have  been  reluctant  to  report  ambiguous 

! results  as  inertial  motions,  and  consequently  their  findings  may  not  be  part  of 

the  literature  on  the  subject. 

According  to  Webster  (1968),  important  early  work  on  the  subject  was  done  by 
Ekman  and  Helland-Hansen,  who  in  1930  made  a series  of  current  measurements  near 
latitude  30°  north.  They  found  diurnal  fluctuations  in  the  signal  that  were 
interpreted  as  tidal  motion.  Ekman  later  concluded  that  these  were  in  fact  inertial 
motions . 

In  later  studies,  Stommel  (1954)  reported  inertial  motions  near  Bermuda,  and 
Swallow  (1957)  found  this  phenomenon  in  the  eastern  Atlantic. 

Investigations  near  latitude  30°  north  during  the  early  1960s  seemed  to 
indicate  that  inertial-period  motions  were  the  result  of  a resonance  phenomenon  of 
the  diurnal  tide  at  the  diurnal  inertial  latitude.  Knauss  (1962)  tracked 
Swallow  floats  at  several  depths  in  the  North  Pacific  near  latitude  28°  north  and 
noticed  a diurnal  clockwise  rotary  motion.  Because  the  recorded  speeds  were 
several  times  greater  than  the  expected  tidal  currents  in  the  open  ocean,  he 
reasoned  that  they  were  due  to  internal  waves  of  inertial  period.  Knauss  con- 
sidered the  concept  of  the  earth's  latitude  bands  acting  as  tuned  circuits,  each 
tuned  to  its  inertial  period.  In  the  vicinity  of  latitudes  29°  and  75°,  the  ocean 
I is  tuned  to  the  diurnal  and  semidiurnal  components,  respectively,  of  the  tide- 

producing  forces.  Strong  inertial  motions,  then,  should  be  found  only  at  those 
latitudes  in  tune  with  generating  forces,  and  nowhere  else. 

A study  by  Reid  (1962)  at  latitude  30°  north  in  the  Pacific  added  support  to 
Knauss'  findings.  Reid  found  that  surface  currents  having  a period  of  rotation 
of  24  hours  were  in  phase  with  the  diurnal  tide.  He  further  noted  an  oscillation 
of  the  temperature  structure  of  12.4-hour  period;  he  suggested  that  the  oscilla- 
tion might  represent  internal  waves  of  lunar  semidiurnal  period. 

Observations  by  Hunkins  (1967)  of  inertial  oscillations  of  ice  island  T-3  in 
the  Arctic  Ocean  tended  to  repudiate  Knauss'  and  Reid's  findings.  Hunkins  noted 
that  "inertial  oscillations  of  floating  ice  are  generated  in  response  to  changing 
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Figure  1.  Study  region  and  locations  of  instruments.  Contours  show  depth  in  meters. 


wind  stress  and  would  presumably  involve  only  the  Ice  and  the  upper  layers  of  the 
ocean.  Their  amplitude  would  reflect  clianges  in  the  wind,  and  their  phase  would  be 
expected  to  change  with  each  wind  system."  His  observed  oscillations  were  of  a 
shorter  period  than  the  semidiurnal  tidal  period,  changed  amplitude  in  response  to 
the  wind,  arbitrarily  changed  phase,  and  were  restricted  to  the  ice  and,  possibly, 

• he  uppermost  water  layers. 

A long  record  of  currents  was  obtained  by  Bernstein  (1972)  from  current  meters 
set  below  ice  island  T-3  in  the  Arctic  Ocean.  He  identified  inertial  motions  in 
the  mixed  layer  below  the  ice  despite  contamination  of  the  current  record  by  the 
semidiurnal  tidal  period,  which  was  nearly  the  same  as  the  local  Inertial  period. 

The  mechanism  of  generation  in  this  case  was  the  stress  produced  by  the  ice  as  it 
moved  over  the  well-mixed  surface  layer.  Among  the  evidence  he  cited  was  the  fact 
that  the  currents  were  much  too  strong  to  be  caused  solely  by  the  tide  and  that 
their  spectral  peak  was  centered  on  a ft. quency  closer  to  the  inertial  than  to  the 
semidiurnal  tidal  frequency. 

Among  the  most  recent  work  on  inertial  motions  is  that  of  Pollard  (1970)  and 
Pollard  and  Millard  (1970),  who  were  able  to  model  inertial  oscillations  of  near- 
surface currents  induced  by  wind  stress  alone.  They  were  successful  in  reproducing 
observed  oscillatory  motions  with  a model  that  assumes  a surface  wind  stress  instan- 
taneously and  uniformly  distributed  as  a body  force  through  the  mixed  layer. 

Perkins  (1970)  observed  inertial  oscillations  in  the  Mediterranean.  He  sepa- 
rated the  inertial  motions  from  the  current  recc  rd  and  determined  their  rate  of 
decay  by  the  method  of  complex  demodulation  to  ju  described  later  in  this  paper. 


II.  ANALYSIS  OF  CURRENTS 


A.  Progressive  Vector  Diagram 

When  the  February  1974  observed  currents  are  plotted  as  a progressive  vector 
diagram,  some  interesting  features  of  the  current  field  emerge.  The  progressive 
vector  diagram  is  simply  the  observed  current  vectors  plotted  end  to  end.  It 
describes  the  path  taken  by  a particle  under  the  influence  of  a horizontally  uni- 
form current  field  identical  to  the  currents  at  the  study  site.  In  examining  the 
progressive  vector  diagram,  it  should  be  remembered  that  it  shows  displacements 
caused  only  by  currents  at  the  measuring  site.  If,  for  example,  a drogue  were 
placed  at  the  study  site  and  tracked  for  a period  of  time,  it  would  encounter  a 
variety  of  current  regimes  and  its  displacement  would  reflect  particle  displacement 
under  the  influence  of  those  regimes.  In  general,  its  path  would  be  different  from 
that  shown  in  Figure  2. 

Figure  2 suggests  rotational  motion  superposed  on  translational  motion.  Espe- 
cially interesting  are  the  loops  in  the  plot  for  the  approximate  periods  4 to  7, 

9 to  15,  17  to  20,  and  25  to  27  February.  This  is  precisely  the  motion  to  be 
expected  when  the  observed  current  has  a significant  rotary  component. 

Examination  of  the  wind  record  at  the  Chevron  rig  (Fig.  3)  indicates  that  the 
observed  rotation  is  associated  with  significant  wind  events.  It  will  be  shown  in 
this  paper  that  the  periods  indicated  above  are  examples  of  inertial  oscillations 
of  the  surface  layer  of  the  ocean.  Other  examples  also  exist  in  the  record  but  are 
not  as  obvious  in  the  progressive  vector  diagram  as  the  above. 
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Figure  3.  Wind  record  at  Chevron  rig  for  4-28  P'ebruary  1974.  F indicates  approxi- 
mate times  of  frontal  passage. 


B.  Data  Filtering 


Extracting  tidal  currents  from  a current  record  can  present  some  annoying 
difficulties.  The  problem  is  one  of  selecting  a filter  or  set  of  filters  that  will 
give  the  truest  possible  reproduction  of  a signal  at  the  desired  frequencies  while 
suppressing  the  signal  at  unwanted  frequencies,  i.e.,  selecting  filters  having  a 
desirable  frequency  response.  The  frequency  response  of  a filter  is  the  ratio  of 
the  amplitude  of  a wave  of  a given  frequency  in  the  time  series  after  filtering  to 
the  original  amplitude.  It  can  be  determined  by  examining  the  effect  of  the  filter 
on  a unit-amplitude  sinusoidal  wave  of  frequency  o (Holloway,  1958).  Such  a wave 
may  be  visualized  as  a unit  vector  in  the  complex  plane  rotating  at  angular  velocity 
2no.  After  application  of  the  filter,  the  modified  wave  may  be  represented  by  a 
new  vector  having  the  same  angular  velocity  as  the  original  but  different  magnitude 
and  phase  angle.  The  frequency  response  of  the  filter  is  then  simply  the  ratio  of 
the  magnitude  of  the  modified  vector  to  t'ii._  unit  vector. 

The  filter  or,  more  precisely,  the  filtering  function  is  a mathematical  opera- 
tor that  consists  of  a series  of  fractional  values  called  weight.^  These  weights 
determine  the  relative  contribution  of  each  observation  in  the  time  series  to  the 
estimate  of  the  filtered  value.  The  weights  are  convolved  with  the  time  series  to 
perform  the  filtering  operation.  For  example,  the  Hamming  filtering  function,  con- 
sisting of  weights  1/4,  1/2,  1/4,  is  symmetric  and  replaces  a time  series  observa- 
tion with  an  average  obtained  from  half  its  value  plus  one-fourth  of  the  two  adja- 
cent values. 

The  frequency  response  of  a filtering  function,  w(t) , is  given  by  the  magni- 
tude of  its  inverse  Fourier  transform,  R(o),  which  is  defined  as 


R(o)  = / w(t)  exp  i(2not)  dt 
^00 

OO  00 

= / w(t)  cos(2Ilot)  dt  + if  w(t)  sin  (2IIat)  dt 


(8) 


where  a is  frequency,  t is  time,  and  i = /-I.  Taking  the  absolute  value, 
I , ( 9 2 I 1/2 

I R(a)  1 = -JlReCR  + [Im(R  (o))]^  • 


(9) 


The  phase  shift,  j),  at  a particular  frequency  produced  by  the  filtering  process 
is  given  by 


c|>  = arc  tan  [ lm(R(a) ) /Re  (R(o) ) ] . 


(10) 


If  the  filtering  function  is  composed  of  discrete  weights,  then  equation  8 
becomes 


m m 

R(o)  =1  Wj^  cos(2IIak)  + i I Wj^  sin(2IIak)  (11) 

k+-n  k=-n 

where  a is  in  cycles  per  data  interval. 
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An  important  requirement  in  the  choice  of  a filter  is  that  it  does  not  shift 
the  phase  of  the  record  at  any  particular  frequency.  This  is  accomplished  by 
choosing  a set  of  weights  that  are  symmetrical  about  some  central  weight,  w^,  i.e., 
w_j  = w. , where  j is  an  integer.  In  this  case,  the  imaginary  term,  the  sine  term 
in  equaiions  8 and  11,  is  zero,  and,  from  equation  10,  41  = 0.  The  frequency 
response  for  the  discrete  case  is  then 

n n 

R(o)  = E Wj^  cos(2nok)  = w^  + 2 E Wj^  cos  (2nok)  (12) 

k=-n  k=l 


for  the  filtering  function  having  (2n  + 1)  weights.  Equation  12  describes  the 
response  of  a filter  with  zero  phase. 

A high-pass  filter  may  be  constructed  very  simply  by  using  a low-pass  filter. 
If  a time  series  is  first  low-passed  and  the  resulting  values  are  then  subtracted 
from  the  original  series,  the  new  series  produced  is  a high-passed  version  of  the 
original  series,  that  is,  one  containing  high-frequency  components  only.  The  fre- 
quency response  of  this  process  may  be  represented  by 


— 1 — (13) 

where  the  subscripts  0 and  LP  denote  the  original  and  low-passed  series,  respec- 
tively. 

Similarly,  a band  pass  of  a time  series  may  be  obtained  by  applying  two  filters 
with  different  frequency  response  to  the  original  data  and  then  subtracting  one 
filtered  series  from  the  other.  The  two  filters  are  chosen  according  to  which  fre- 
quencies are  to  be  suppressed.  Application  of  one  filter  is  used  to  obtain  a 
record  with  only  low  and  intermediate  frequencies.  The  second  filter  is  used  to 
pass  only  low  frequencies.  The  result  of  subtracting  the  second  series  from  the 
first  is  a band  pass,  which  may  be  expressed  in  terms  of  the  frequency  response  of 
the  filters  as 


R(o)lpi-lP2  ~ ^^'^^LPl  (1^) 

where  LPl  and  LP2  represent  low-pass  filters  1 and  2,  respectively.  This  process 
will  remove  both  high-  and  low-frequency  components  and  will  yield  values  at  inter- 
mediate frequencies. 

The  problem  of  extracting  tidal  currents  from  a current  time  series  is  one  of 
appropriately  band-pass  filtering  that  series.  An  effective  and  surprisingly  sim- 
ple filter  for  analysis  of  tides,  which  was  used  in  this  study,  was  developed  by 
Doodson  and  Warburg  (1941).  The  weights  of  the  Doodson-Warburg  filter,  which  are 
used  to  multiply  a time  series  with  At  = 1 hour,  are  shown  in  Table  1.  These 
weights,  which  Doodson  and  Warburg  call  "multipliers  for  mean  sea  level,"  were 
derived  empirically  by  combining  several  simple  filters  that  discriminate  against 
the  principal  tidal  constituents  and  higher  frequencies  (Doodson  and  Warburg,  1941; 
Groves,  1955).  The  Doodson-Warburg  filter  has  the  desirable  property  of  being 
symmetric  and  therefore  does  not  change  the  phase  of  the  time  series. 
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Table  1 

Numerators  of  Weights  of  the  Doodson-Warburg  Filter 


(1) 

1 

(9) 

1 

(17) 

1 

(25)  0 

(33) 

0 

(2) 

0 

(10) 

0 

(18) 

1 

(26)  1 

(34) 

1 

(3) 

1 

(11) 

2 

(19) 

2 

(27)  1 

(35) 

0 

W 

0 

(12) 

0 

(20) 

0 

(28)  0 

(36) 

0 

(5) 

0 

(13) 

1 

(21) 

2 

(29)  2 

(37) 

1 

(6) 

1 

(14) 

1 

(22) 

1 

(30)  0 

(38) 

0 

(7) 

0 

(15) 

0 

(23) 

1 

(31)  1 

(39) 

1 

(8) 

1 

(16) 

2 

(24) 

2 

(32)  1 

NOTE: 

Each  value 

is 

divided  by  30  to  i 

obtain  the 

set 

of  weights. 

Another 

filter  used 

in 

this  analysis  is 

the  31-point 

binomial  filter. 

It  also 

is 

a symmetrical  filter 

whose  weights  are  proportional 

to 

the  coefficients 

of  the 

expansion  of 

the  binomial  form  (p  + q)*^,  where  n is  the  number  of  desired 

weights 

less  1 (n  = 

30  in  this  case 

) (Holloway,  1958) 

. In  general. 

the  frequency  response 

of 

a binomial  filter  as 

calculated  from  equation  8 is 

given  by 

R(o)  = cos’’  (noAt) 


(15) 


where  At  = 1/3  hour  in  the  present  case. 

Recall  that  for  increasing  n,  the  envelope  of  the  coefficients  of  the  expan- 
sion of  (p  + q)"  approaches  a Gaussian  curve.  In  both  filters  the  sum  of  the 
weights  is  unity.  This  property  of  the  two  filters  is  particularly  important  inas- 
much as  each  weight  expresses  in  what  proportion  each  of  a set  of  observations  is 
to  be  counted.  Each  weight  is,  therefore,  a fractional  part  of  a sum  that  is  1. 

The  original  current  record  was  decomposed  into  vector  components  and  filtered 
using  the  binomial  filter  in  order  to  remove  the  high  frequencies.  The  original 
values  were  also  low  passed  with  the  Doodson-Warburg  filter  to  remove  the  principal 
tidal  constituents.  Because  the  current  values  were  available  at  intervals  of 
1/3  hour  and  the  Doodson-Warburg  filter  requires  a 1-hour  interval,  the  filter  was 
first  applied  at  1-hour  intervals  beginning  at  t^  = 0 hour.  It  was  then  time 
shifted  by  1/3  hour  and  applied  at  1-hour  intervals  beginning  with  starting  time 
tp  = 1/3  hour  and  again  with  t^  = 2/3  hour.  The  result  is  a filtered  time  series 
with  At  = 1/3  hour.  The  output  of  the  Doodson-Warburg  filter  was  then  subtracted 
from  the  output  of  the  binomial  filter  to  obtain  the  vector  components  of  the 
tidal  currents.  In  terms  of  the  frequency  response,  the  process  may  be  expressed 
as 


’^(®)bD-PASS  “ " ’^(o)dw 

where  BD-PASS  represents  the  band-pass  process,  BN  is  the  binomial  filter,  and  DW 
is  the  Doodson-Warburg  filter. 
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Figure  A is  a plot  of  the  frequency  response  functions  of  the  filters  dis- 
cussed. Notice  that  the  31-point  binomial  filter  eliminates  all  but  about  10  per- 
cent of  the  energy  at  0.35  CPU,  yet  it  preserves  about  90  percent  of  the  semidiurnal 
tide  and  nearly  all  of  the  diurnal  tide. 

On  the  other  hand,  the  Doodson-Warburg  filter  is  highly  effective  in  eliminating 
the  lower  frequencies,  specifically  the  diurnal  and  semidiurnal.  Table  2,  adapted 
from  Groves  (1955),  gives  the  frequency  response  of  the  filter  for  the  principal 
tidal  constituents.  The  high-pass  filter  formed  from  the  Doodson-Warburg  filter, 
therefore,  preserves  almost  completely  those  motions  associated  with  tidal  forcing. 
The  band-pass  signal  found  using  the  binomial  and  Doodson-Warburg  filters  contains 
about  97  percent  of  the  current  magnitude  near  the  diurnal  period,  90  percent  near 
the  semidiurnal,  and  65  percent  near  the  quarter-diurnal.  Energy  at  periods  greater 
than  2 days  is  greatly  reduced. 

A band-passed  version  of  the  February  data  was  obtained  and  plotted  in  Figure  5 
as  hodographs  for  each  2A-hour  period  beginning  at  0000  hour  and  ending  at  2300  hour 
of  each  day.  The  hodograph  represents  the  path  traced  by  the  tip  of  the  tidal  cur- 
rent vector  over  a 24-hour  period.  If  we  assume  that  large  tidal  range  (Fig.  6) 
must  produce  strong  tidal  currents  and  small  tidal  range  must  produce  weak  tidal 
currents  (Sverdrup  et  al.,  1942),  the  relative  sizes  of  the  hodographs  are  rather 
surprising.  For  example,  31  January  through  3 February  is  a period  of  rather  large 
tidal  range,  but  the  hodographs  indicate  generally  weak  currents.  From  8 through  12 
February  is  a period  of  equatorial  tides  (minimum  tidal  range),  but  the  hodographs 
are  the  largest  of  the  study.  Consider  also  the  hodographs  of  25  and  26  February, 
which  show  strong  tidal  currents  for  a period  of  equatorial  tides. 


Table  2 

Frequency  Response  of  Doodson-Warburg  Filter 
to  the  Principal  Tidal  Constituents 


Tidal  Constituent  Period  R(a) 

(Hours) 


M2 

12.42 

-0.00058 

4 

12.00 

0.00000 

N2 

12.66 

0.00171 

K2 

11.97 

0.00033 

Kl 

23.93 

0.00015 

0l 

25.82 

0.00299 

Pi 

24.07 

-0.00013 

It  appears  that  some  process  other  than  astronomical  tidal  forcing  is  putting 
energy  into  the  tidal  frequency  band.  The  wind  record  (Fig.  3)  obtained  at  the 
Chevron  platform  and  the  Weather  Bureau  record  for  Boothville,  Louisiana  (not  shown) 
reveal  a very  strong  correlation  between  times  of  shifting  wind  direction  and 
increased  wind  speed  (usually  associated  with  atmospheric  fronts)  and  large  well- 
defined  "tidal  ellipses."  At  other  times  the  hodographs  have  a generally  ragged 
and  somewhat  confused  appearance. 
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31 -pt  BN  Filter 
D-W  Filter 


Figure  4.  Frequency  response  of  binomial  filter  (BN),  Doodson— Warburg  filter  (DW) , and  band— pass 
process  (BD-Pass).  D,  f,  and  SD  represent  the  diurnal,  inertial,  and  semidiurnal  frequencies, 
respectively. 


Figure  6.  A.  Predi..ted  tidal  elevations  at  Caminada  Pass  from  National  Ocean  Survey  tables,  Department  of  Commerce. 
B.  Measured  tidal  elevations  at  Bayou  Rigaud  collected  by  U.S.  Coast  Guard,  Grand  Isle.  (Note  different  scales  on 
ordinate  axes.)  Symbols  along  horizontal  axis  are  as  follows:  e = first  lunar  quarter;  N = moon  farthest  north  of 
equator;  P = moon  at  perigee;  o = full  moon;  E = moon  on  equator;  a = last  lunar  quarter;  S = moon  farthest  south  of 
equator;  A = moon  at  apogee;  • = new  moon.  Dates  are  for  January  and  February  197A. 


The  times  of  frontal  passage  were  determined  using  National  Weather  Service 
weather  maps;  wind  shifting  at  the  Chevron  rig,  Boothv'lle,  and  New  Orleans;  dry 
bulb  and  wet  bulb  temperatures  at  Boothvllle;  atmospheric  pressure  at  Boothville; 
and  times  of  precipitation  at  Boothville  and  New  Orleans. 

Nearly  every  wind  shift  produces  large  hodographs,  which  diminish  in  size 
through  a period  of  1 to  4 days.  A front  on  3 February  was  followed  by  large  hodo- 
graphs on  4 through  6 February.  A wind  shift  on  9 February  was  followed  by  several 
days  of  large  shrinking  hodographs  during  equatorial  tides.  On  18  February  a front 
produced  fairly  well  behaved  rotary  currents  for  a day.  The  wind  shift  of  20 
February  produced  a nearly  circular  hodograph.  The  front  on  21  February  and  wind 
shift  on  the  22nd  produced  a well-defined  hodograph.  Finally,  a front  on  24  Febru- 
ary seems  responsible  for  the  strong  rotary  currents  that  follow  on  25  and  26 
February,  a time  of  equatorial  tides.  Possibly  the  only  exception  to  this  pattern 
is  the  hodograph  of  16  February,  which  follows  the  front  of  15  February  and  shows 
rather  poorly  behaved  clockwise  rotary  motion. 

The  evidence  suggests  that  band-passed  data  contain  rotary  currents  that  are 
induced  partially,  if  not  in  the  most  part,  by  wind.  In  the  reminder  of  this  paper 
an  attempt  will  be  made  to  demonstrate  that  these  are  indeed  inertial  oscillations 
of  the  surface  layer. 


C.  Complex  Demodulation  of  Observed  Currents 

A rotating  current  having  both  clockwise  and  counterclockwise  components  at 
many  different  frequencies  may  be  expressed  as 


00 

u(t)  = / A(uj)  exp  (-iut)  do) 


(17) 


where  A(u))  is  the  complex  amplitude  of  the  current  at  radian  frequency  m and  t is 
time.  For  the  case  of  a current  signal  dominated  by  inertial  oscillations,  we 
assume  that  A(u))  is  sharply  peaked  near  to  = f . Complex  demodulation  (Tukey,  1961; 
Perkins,  1970;  Bloomfield,  1976)  of  this  signal  using  f,  the  theoretical  inertial 
frequency,  as  the  demodulating  frequency  yields 

U(t)  = u(t)  exp  ( ift) 

00 

= / A(oj)  exp  (-imt)  exp  ( 1ft)  d^ 

—GO 

CO 

= I A(u)  exp  [ i(f  - u))t  ] d|^  . (18) 


In  trigonometric  form,  equation  18  is 

00 


U(t)  = 

Low-pass  filtering  of 


/ A(oj)  [cos  (f  - (jj)t  + i sin(f  - 

— CD 

this  signal  will  have  one  of  two 


iv)t]  djj  . 

effects : 


(19) 
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1.  If  u)  is  close  to  f,  the  amplitude  of  the  signal,  A(a>) , will  be  relatively 
unchanged  because  the  signal's  frequency  will  be  low. 

2.  If  u)  is  not  close  to  f,  a high-frequency  signal  will  result  and  it  will 
therefore  be  eliminated. 

The  result  of  low-pass  filtering  yields  the  demodulated  signal 

CD 

Ujj(t)  = / A(m)  G(e)  exp  ( iet)  d^  (20) 


where  G(e)  is  the  transfer  function  of  the  low-pass  filter  and  c = f - w.  It  is 
desirable  to  use  a filter  that  will  not  alter  the  phase  of  the  output  signal.  A 
two-sided  filter  having  symmetrical  weights  has  such  a property.  One  such  filtering 
function  is  the  equally  weighted  "running  mean"  or  "moving  average,"  whose  transfer 
function  (Holloway,  1958)  is  given  by 


G(m)  = (sin  ujT)/wT  (21) 


where  T is  the  half-width  of  the  filter  in  the  time  domain.  Thus,  for  <ij->  f,  £ ->•  0, 
G(e)  ■*  1,  exp  (iet)  1,  and  0]j(t)  -►  A(f). 

What  the  technique  accomplishes,  then,  is  to  transform  a signal  at  frequency  w, 
which  may  be  concealed  by  large  contributions  at  other  frequencies,  to  a demodulated 
signal  having  a much  reduced  frequency  of  e,  which,  when  plotted,  is  more  readily 
apparent  to  the  eye.  The  demodulated  current  has  a phase  implicit  in  the  complex 
amplitude  A(u)  defined  as  if  = tan"i  \Im  [A(m)]/Re  [A(m)]}. 

Complex  demodulation  may  also  be  thought  of  as  a local  version  of  harmonic 
analysis  in  that  it  is  a means  to  describe  a signal  at  a given  frequency  in  terms  of 
its  amplitude  and  phase.  It  is  local  in  the  sense  that  it  describes  these  param- 
eters in  the  neighborhood  of  a particular  time,  t,  rather  than  for  the  entire 
record.  In  this  study,  a filtering  interval  of  two  inertial  periods  was  chosen, 
i.e.,  T = 1 inertial  period.  Therefore,  each  Un(t)  computed  is  a function  only  of 
those  currents  observed  over  the  time  period  t f 1 inertial  period  and  is  indepen- 
dent of  all  other  observed  current  velocities. 

Demodulation  of  the  February  current  record  yields  some  Interesting  results 
(Fig.  7)  when  compared  with  the  wind  record.  At  first  glance,  there  is  an  excellent 
correlation  between  rapid  shifts  in  the  wind  direction  (usually  associated  with 
fronts)  and  peaks  in  the  amplitude  of  the  demodulated  current.  Furthermore,  these 
peaks  are  generally  accompanied  by  a flattening  of  the  plot  of  the  phase  of  the 
demodulated  current,  indicating  a signal  of  inertial  or  near-inertial  frequency. 

Because  the  inertial  and  diurnal  tidal  periods  are  so  close,  tidal  currents 
will  unfortunately  also  be  preserved  in  the  demodulated  signal  with  little  change 
in  magnitude.  However,  because  the  peaks  in  current  magnitude  on  Figure  7 all 
appear  to  be  associated  with  specific  wind  events,  it  seems  unlikely  that  tidal 
currents  contribute  substantially  to  the  demodulated  current.  Furthermore,  if  the 
peaks  in  magnitude  were  tidal  rather  than  inertial,  then  the  associated  phase 
would  have  a constant  rate  of  change  associated  with  each  magnitude  peak  and  pre- 
sumably for  the  entire  record.  This  is  clearly  not  the  case  in  Figure  7 or  8. 
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Figure  7.  February  demodulated  current.  The  F indicates  approximate  times  of  frontal  passage 


Thesse  points  will  be  discussed  further  in  connection  with  specific  peaks  in  the 
current  magnitude. 

There  are  tliree  prominent  broad  peaks  in  tlie  amplitude,  centered  on  5,  11,  and 
26  February,  that  have  remarkably  constant  pliase.  Some  of  the  less  prominent  peaks 
have  a phase  that  is  slowly  changing;  a rotation  of  the  current  vector  either  faster 
or  slower  than  inertial  is  suggested.  An  increasing  phase,  as  defined  tiere,  indi- 
cates that  the  current  vector's  period  of  rotation  is  less  than  inertial,  whereas 
a decreasing  phase  results  from  a vector's  rotating  with  a period  greater  than 
inertial . 

There  are  two  possible  features  in  a wind  field  that  are  very  effective  in 
generating  or  changing  the  amplitude  of  inertial  oscillations:  (1)  a strong  wind 

accompanied  by  a fairly  sudden  shift  in  direction  or  (2)  a wind  vector  rotating 
clockwise  for  up  to  half  an  inertial  period  (Pollard  and  Millard,  1970;  Pollard, 
1970).  As  an  example  of  case  1,  consider  a wind  stress  acting  in  one  direction  for 
an  arbitrary  period  of  time.  Such  a wind  adds  momentum  to  an  inertial  current 
during  the  first  half  of  each  inertial  period  and  removes  it  during  the  second  half. 
Therefore,  a wind  that  suddenly  changes  direction  or  is  removed  at  the  end  of  1/2, 
3/2,  5/2,  etc.,  inertial  periods  produces  maximum  inertial  amplitudes.  For  case  2, 
if  a wind  vector  rotates  clockwise  at  inertial  frequency,  it  continuously  adds 
momentum  to  an  inertial  current,  thereby  continuously  increasing  its  amplitude. 

Both  conditions  are  usually  present  in  varying  degrees  when  an  atmospheric 
front  approaches  or  passes  a given  site.  Along  the  Louisiana  coast,  an  approaching 
front  is  normally  accompanied  by  a relatively  slow  clockwise  rotation  of  the  wind, 
which  blows  first  from  the  east,  then  from  the  south,  and,  finally,  as  the  front 
passes,  suddenly  shifts  to  the  west  or  northwest.  The  passage  of  the  front  satis- 
fies case  1,  whereas  wind  conditions  accompanying  the  approach  of  the  front  satisfy 
case  2. 

Another  example  of  case  2 can  occur  without  the  presence  of  a front  as  a result 
of  the  geometry  of  the  isobars.  For  example,  a high-pressure  system  whose  center 
passes  north  of  an  oceanic  site  toward  the  east  has  clockwise  rotating  winds. 
Depending  on  the  rate  at  which  the  system  moves,  the  rotation  rate  of  the  accompany- 
ing winds  may  approach  inertial  frequency.  However,  a clear  example  of  this  situa- 
tion is  not  present  for  the  time  period  of  this  study. 

Another  factor  that  is  crucial  in  determining  the  amplitude  of  the  wind-induced 
oscillations  near  30°  latitude  is  the  magnitude  and  direction  of  the  diurnal  tidal 
current.  If,  for  example,  the  tidal  current  has  a component  opposite  a wind  vector 
component  for  any  period  of  time,  momentum  will  be  removed  from  the  current  field. 
The  result  will  be  diminishing  or  destruction  of  rotary  motion.  If,  on  the  other 
hand,  the  wind  and  tidal  current  vectors  each  have  a component  acting  in  the  same 
direction,  momentum  will  be  added  to  the  current  and  rotary  motion  will  be  enhanced. 
The  same  arguments  hold  for  pre-existing  inertial  motion  at  the  time  of  a wind 
event . 

It  is  instructive  to  compare  wind  strength  and  duration  with  the  amplitude  of 
the  current  peaks.  For  example,  the  peaks  of  11  and  18  February  are  associated 
with  generating  winds  of  approximately  equal  strength  but  different  duration.  The 
wind  on  9 February  peaks  sharply  and  rotates  at  inertial  frequency,  and  truly  strong 
winds  occur  for  about  12  hours.  The  wind  of  18  February  shifts  rapidly  but  becomes 
unidirectional  and  still  strong  on  19  February,  and  the  peak  lasts  about  24  hours. 
Results  of  the  two  events  are  quite  different.  The  wind  event  of  9 February 
increases  the  already  existing  oscillations  by  about  20  cm/sec,  whereas  the  peak 
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produced  on  18  Kebrnary  is  only  7 cm/sec  above  the  background  and  is  eliminated  by 
tbe  strong  unidirectional  wind  of  19  February. 

Winds  following  tbe  front  of  24  February  produce  some  interesting  amplification 
as  well  as  destruction  of  inertial  motion.  Following  the  wind  shift  in  the  early 
morning  of  24  February,  inertial  current  amplitude  built  for  half  the  inertial  per- 
iod and  then  diminished  for  tlie  second  half  of  the  period.  Very  slow  rotation  of  a 
strong  wind  vector  following  the  wind  shift  seems  to  be  responsible.  Strong  ampli- 
fication of  tile  inertial  current  during  25  February  and  early  26  February  for  some- 
wliat  less  tlian  an  inertial  period  was  caused  by  a moderate  to  strong  wind  vector  that 
rotated  slowly  clockwise.  Because  February  is  a time  of  equatorial  tides,  lunar 
effects  are  minimal  and  tliercfore  the  motion  must  be  considered  inertial.  There 
is  some  destruction  and  tlien  amplification  for  about  1/2  inertial  period  following 
the  maximum  current  of  early  26  February.  The  ptiase  of  the  demodulated  current  shows 
a nearly  constant  value  from  early  25  February  to  the  end  of  the  record;  a frequency 
very  close  to  inertial  is  indicated. 

According  to  the  meteorological  data  from  the  Weather  Bureau  station  at  Booth- 
ville,  the  front  of  3 February  is  accompanied  by  a strong  narrow  peak  in  the  wind 
speed  and  followed  by  several  days  of  relative  calm.  (Unfortunately,  the  anemometer 
at  the  Chevron  rig  was  not  activated  until  4 February).  There  is  an  associated 
strong  uniform  increase  in  the  demodulated  current  with  a constant  phase.  Tidal 
effects  might  seem  significant  because  of  closeness  to  the  time  of  tropic  tides. 
However,  the  moon  Is  at  its  maximum  declination  on  2 February  and  the  demodulated 
current  for  that  day  is  only  10  cm/sec  as  compared  with  the  peak  of  40  cm/sec  on 
5 February,  following  the  front  of  3 February.  The  lime  lag  between  maximum  lunar 
declination  and  resulting  tropic  tide,  the  so-called  "age"  of  the  tropic  tide,  along 
the  Louisiana  coast  is  about  6 hours  (Marmer,  1954).  As  a result,  the  peak  of 
5 February  occurs  much  too  late  to  be  the  result  of  the  lunar  effects  of  2 February. 

Peaks  in  the  demodulated  current  for  the  period  15  to  24  February  are  relatively 
small  and  of  short  duration.  From  the  nature  of  the  wind  during  this  time  interval, 
it  appears  that  wind  stress  is  as  important  a factor  in  the  decay  of  the  currents  as 
it  is  in  their  generation.  Each  peak  has  a decay  time  of  12  to  24  hours  which  in 
each  case  corresponds  to  a period  of  fairly  unidirectional  wind  lasting  12  to  24 
hours.  Consider,  for  example,  the  peak  centered  on  23  February.  The  wind  shift 
that  begins  on  22  February  is  followed  on  23  February  by  light  winds  toward  the 
west.  The  inertial  current  on  23  February  shows  a uniform  decrease  in  amplitude. 

This  decrease  suggests  chat  wind  stress  is  tending  to  remove  momentum  from  the 
current  system. 

Furthermore,  the  close  spacing  of  strong  wind  episodes  during  this  period  seems 
to  preclude  the  observation  of  long-duration,  strong  inertial  amplitudes  in  the 
current.  Spacing  of  wind  shifts  and  fronts,  with  the  exception  of  the  one  on  15 
February,  is  about  one  to  two  inertial  periods,  a time  interval  that  is  rather  short 
when  one  wislies  to  resolve  well  this  type  of  rotary  motion. 

In  order  to  further  substantiate  the  conclusions  drawn  from  the  demodulation 
of  the  February  current  record,  a current  time  series  extending  from  25  .April 
through  2 .June  was  tIso  subjected  to  complex  demodulation.  Figure  8 shows  only  one 
major  peak  in  the  demodulated  current  centered  on  12  May  and  associated  with  a 
frontal  passage  early  on  11  May  that  exhibited  gusts  of  up  to  21  m/sec.  Although 
weather  maps  for  this  time  period  indicate  the  passing  of  four  fronts,  weather  data 
from  Boothville  indicate  that  the  front  of  11  May  was  by  far  the  strongest,  the 
others  being  fairly  diffuse  and  generally  barely  perceptible  in  the  study  region. 
There  is  some  correspondence  between  the  shifting  of  the  wind  at  Boothville  and 
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tile  small  current  peaks,  but  it  is  difficult  to  characterize  these  as  inertial  owing 
to  their  relative  weakness.  in  fact,  this  record  seems  to  be  indicative  of  back- 
ground tidal  noise  for  this  region.  On  the  basis  of  this  record,  a fairly  repre- 
sentative value  for  maximum  tidal  current  magnitude  is  on  the  order  of  10  cm/sec. 


D.  Spectral  Anal^’^1^ 

The  Fourier  transform  (Blackman  and  Tukey,  1959)  of  a finite-duration  time 
series  x(t)  is  given  by 


T/2 

X((d)  = - / x(t)  e"i“t  dt,  -T/2  - t - T/2  (22) 

T _t/2 

where  T is  the  duration  of  the  series,  o is  radian  frequency,  and  t is  time. 

The  correlation  function  for  components  u and  v of  a finite-duration  vector 
time  series  is  given  by 


T/2 

= - / u(t)v(t  + t)  dt 
T _t/2 


(23) 


where  t denotes  a time  lag.  The  symbol  represents  the  cross  correlation.  Simi- 
larly, Rjjy  and  R^y  may  be  used  to  express  autocorrelations  of  u and  v,  respectively. 

The  two-sided  spectrum  is  defined  by  the  Fourier  transform  of  the  correlation 
function 


S (u) 
uv  ' 


/ R„^(t)  e-i“"  dx 


(24) 


for  -o>  < m < “>.  This  quantity  is  referred  to  as  the  cross-spectrum  of  series  u and 
V.  We  may  similarly  compute  the  autospectrum  of  u,  from  Ryu,  and  Syvi  the 

autospectrum  of  v,  from  Rw-  These  are  commonly  termed  the  power  spectra  of  u and 
V,  respectively. 

Conventional  analysis  of  time  series  usually  deals  with  the  so-called  "physi- 
cally realizable  one-sided  spectra,"  which,  as  defined  by  Bendat  and  Piersol  (1966), 
are 


Pw(w)  = 2Svv(>^)  (25) 


^uv(*^^  - 2S^j^(u) 


where  the  spectra  are  now  defined  for  0 < co  < <*>.  The  cross  spectrum  can  further 
be  divided  into  a real  part,  C^^,  called  the  cospectrum,  and  an  imaginary  part, 
Qyy,  called  the  quadrature  spectrum,  expressed  as 
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(26) 


for  0 < u)  < ■»>. 

According  to  Bendat  and  Piersol  (1966),  is  an  odd  function  of  c;,  i.e., 

Quv(“)  = -Quv(“)>  Pvv(“)  ‘>'■0  even  functions,  i.e.,  Pmj(->.))  = Puu(-) 

and  Pvv(““)  ~ Pvv(‘^)*  Using  the  quantities  in  c,'j"ations  25  and  26,  Gonella  (1972) 
writes  for  the  one-sided  counterclockwise  spectrum 

S+  = ^ [Puu('^0  + Pvv'^**'^ 


For  the  one-sided  clockwise  spectrum, 

S_  = “ [Piiu(‘^)  Pvv^“^  “ 2Q„y(uj)] 


(28) 


The  total  spectrum  is 


St(eO  = S_  + S+ 

= 2 [Puu(‘^)  + PvvC'-)! 

We  may  now  write  the  difference  between  the  two  spectra: 


(29) 


S_  - S+  = -Q^^M 


(30) 


The  sign  of  this  quantity  represents  the  polarization  of  energy  in  the  clockwise 
sense  for  positive  sign  and  counterclockwise  for  negative  (Gonella,  1972). 

The  ratio 


Cr  = (S_  - S+)/Sjm)  = -2Quv(‘^)/(Puu(“)  + Pw(‘‘')  J • 


(31) 


called  the  rotary  coefficient  (Mooers,  1973;  Gonella,  1972;  Fofonoff,  196'^.),  gives 
a measure  of  the  partition  of  total  spectral  energy  S^(,)  between  clockwise  and 
counterclockwise  rotating  components.  Its  value  is  zero  for  unidirectional  motion, 

-1  for  pure  counterclockwise  rotary  motion,  and  +1  for  pure  clockwise  motion. 

Spectra  for  the  data  set  were  computed  with  the  Cooley-Tukcy  Fast  Fourier 
Transform  (FFT).  A potentially  troublesome  aspect  encountered  when  using  the  FFT 
stems  from  the  fact  that  the  time  series  being  analyzed  is  finite.  TIk,  FFT  "ass.jmes" 
that  a given  signal  repeats  Itself  infinitely  in  time  (W.  Wiseman,  1976,  personal 
communi  at  Ion) . Therefore,  if  a given  time  series  contains  sinusoidal  components 
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and  the  record  to  be  analyzed  is  such  that  the  beginning  and  the  end  of  the  record 
do  not  correspond  to  the  same  point  on  the  sine  curve,  the  ends  of  the  record  will 
be  interpreted  as  discontinuities.  A strong  high  frequency  flattening  in  the  spec- 
trum may  be  produced.  However,  the  problem  can  be  circumvented  by  tapering  the  ends 
of  the  time  series.  The  February  current  record  was  linearly  detrended  and  tapered 
over  10  percent  of  the  record  at  the  beginning  and  the  end  of  the  time  series  using 
a cosine  bell,  as  suggested  by  Brigham  (1974). 

The  FFT  algorithm  has  certain  advantages  over  the  computing  of  spectral  esti- 
mates with  the  time-lagging  method.  For  example,  the  number  of  spectral  estimates 
that  it  can  compute  is  equal  to  half  the  number  of  data  points,  whereas  one  is 
restricted  in  the  lagging  technique  to  about  10  percent  of  the  number  of  data 
points.  The  FFT,  therefore,  achieves  considerably  better  frequency  resolution.  Its 
iimitlng  frequency,  the  Nyquist  frequency,  fj^  = l/2At,  where  At  is  the  data  interval 
represents  the  highest  frequency  for  which  spectral  information  may  be  obtained.  An 
important  advantage  in  statistical  analysis  is  the  fact  that  the  estimates  computed 
by  the  FFT  are  nearly  statistically  independent.  Finally,  it  is  highly  efficient 
in  terms  of  computer  time  when  compared  with  other  techniques  for  obtaining  the 
same  amounts  of  information. 


In  using  the  FFT,  we  write  equation  24  in  its  alternate  form  (Lathi,  1968) 


S (oj) 
uv 


1 

T 


.T/2 

-T/2 


u(t) 


1 

T 


,T/2 

-T/2 


v(t) 


-iu)t 

e 


dt 


(32) 


where  the  asterisk  denotes  the  complex  conjugate  of  the  Fourier  transform.  We  next 
define  real  and  imaginary  Fourier  coefficients  a(u))  and  b(oj)  for  u(t)  and  c(uj)  and 
d(co)  for  v(t),  which  may  be  used  to  express  the  cross  spectrum  in  terms  of  trigo- 
nometric functions  (Lathi,  1968).  Therefore,  dropping  the  functional  notation  (to) 
from  a (to),  b(io),  c(to),  and  d(to). 


1 

T 


1 

T 


T/2  . 1 * 

/ u(t)  dt  = (a  + ib)/2 

-T/2  J 

T/2 

f v(t)  e dt  = (c  - id)/2 

-T/2 


(33) 


The  two-sided  cross  spectrum  is  then 


Suv(*^)  = + bd  + i(bc  - ad)]/4 

Finally,  the  one-sided  autospectra  are 


Puu(“) 

= 

(to)  = 

(a2 

+ b2)/2 

Pvv(“) 

= 2Sw 

(to)  = 

(c2 

+ b2)/2 

(34) 


(35) 
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and  one-sided  cospectrum  and  quadrature  spectrum  are,  respectively,  the  real  anf" 
imaginary  parts  of 


C^^(u.)  = 2Re  [S^^(m)]  = (ac  + bd)/2 


Q (oj)  = 21m  [S  (w)]  = (be  - ad)/2 
uv  uv 


(36) 


where  0 < m < “. 

The  various  spectral  quantities  were  computed  using  the  FFT,  and,  from  these, 
values  of  Cr  were  calculated  at  all  frequencies.  Table  3 lists  some  values  of  the 
rotary  coefficient,  as  well  as  the  clockwise  and  counterclockwise  spectral  energy 
for  frequencies  near  the  inertial-diurnal  and  semidiurnal  frequencies.  As  expected, 
the  major  clockwise  spectral  peaks  occur  in  the  inertial-diurnal  region.  These 
particular  frequencies  also  generally  have  a Cr  > 0.9  and  Cr  > 0.99  immediately 
above  the  inertial  and  diurnal  frequencies.  The  lack  of  a stronger  spectral  peak 
in  the  semidiurnal  tidal  frequency  is  rather  surprising  in  view  of  the  fact  that 
the  tide  for  the  Louisiana  coast  is  semidiurnal  for  part  of  the  month.  Also  sur- 
prising are  the  values  of  Cr  in  the  vicinity  of  the  semidiurnal  frequency.  The 
values  indicate  nearly  unidirectional  motion  instead  of  the  expected  clockwise 
polarization.  At  least  for  this  record,  it  must  be  concluded  that  the  semidiurnal 
tidal  component  makes  no  significant  contribution  to  rotary  motion  for  the  region 
being  studied. 

Figure  9 shows  the  counterclockwise,  clockwise,  and  total  spectra.  Except  for 
rather  minor  spikes,  there  is  no  significant  energy  in  any  of  the  spectra  in  the 
vicinity  of  the  semidiurnal  tidal  frequency.  The  counterclockwise  spectrum  contains 
significant  energy  only  at  the  lowest  frequencies;  spectral  energy  drops  off  fairly 
uniformly  for  the  higher  frequencies.  As  expected,  energy  in  the  clockwise  and 
total  spectra  is  concentrated  in  a major  peak  near  the  inertial-diurnal  frequencies. 
By  far  the  strongest  peak  in  S_  and  occurs  at  a frequency  of  0.043  CPH,  corres- 
ponding to  a period  of  23.26  hours.  There  are  four  less  energetic  peaks  in  the 
vicinity  associated  with  frequencies  of  0.031,  0.035,  0.049,  and  0.053  CPH  (periods 
of  32.26,  28.57,  20.41,  and  18.87  hours,  respectively). 

Substantial  energy  is  present  in  all  three  spectra  for  very  low  frequencies. 
This  energy  is  apparently  associated  with  a mean  flow  in  the  region. 

The  strongest  peak,  representing  inertial-tidal  frequencies,  occurs  at  a 
somewhat  higher  frequency  than  theory  indicates  for  inertial  currents  and  is  spread 
over  a band  of  frequencies.  Gonella  (1971)  similarly  observed  a shorter  period  of 
inertial  rotation  over  a band  of  frequencies  at  the  Bouee  Laboratoire  in  the 
Mediterranean.  The  phenomenon  results  from  the  rotary  motion  at  various  depths  in 
the  water  layer  tending  toward  tlie  inertial  frequency. 

Let  us  assume  that  at  some  arbitrary  time,  t^  = 0,  a wind  stress  acting  at  the 
surface  is  suddenly  removed  and  the  resulting  vertical  current  distribution  is  given 
by  Uq(z).  The  vertically  integrated  velocity  distribution,  Ro . at  time  to  is  then 

= L (37) 
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Table  3 


Spectral  Estimates  and  Rotary  Coefficient  near  the  Inertial-Diurnal  and  Semidiurnal 
Frequencies  for  the  February  1974  Current  Record 


Frequency 

(CPH) 

S_ 

((cm/sec)2/CPH) 

((cm/sec)^/CFH) 

S 

0.027 

2.173 

0.449 

0.658 

0.028 

1.110 

0.376 

0.494 

0.030 

2.692 

0.588 

0.642 

0.031 

10.635 

0.  383 

0.930 

0.032 

6.041 

0.305 

0.904 

0.034 

2.515 

0.428 

0.  709 

0.035 

10.908 

0.071 

0.987 

0.036 

5.461 

0.450 

0.848 

0.038 

0.437 

0.  718 

-0.243 

0.039 

2.986 

0.777 

0.587 

0.041 

22.092 

0.423 

0.962 

0.042 

27.264 

0.045 

0.997 

0.043 

53.460 

0.207 

0.992 

0.045 

4.921 

0.013 

0.994 

0.046 

3.427 

0.465 

0.761 

0.047 

6.956 

0.274 

0.924 

0.049 

15.422 

0.473 

0.940 

0.050 

7.561 

0.383 

0.  904 

0.051 

6.517 

0.043 

0.987 

0.053 

10.965 

0.139 

0.975 

0.054 

0.064 

0.076 

-0.087 

0.055 

2.298 

0.206 

0.836 

0.076  0.061  0.228  -0.578 


0.077 

0.431 

0.193 

0.  382 

0.078 

0.254 

0.150 

0.258 

0.080 

0.172 

0.205 

-0.087 

0.081 

0.044 

0.044 

-0.005 

0.082 

0.432 

0.414 

0.021 

0.084 

0.126 

0. 135 

-0.032 

0.085 

0.412 

0.183 

0.  385 

0.086 

0.006 

0.004 

0.213 

0.088 

0.129 

0.058 

0.  380 

0.089 

0.078 

0.185 

-0.407 

0.091 

0.099 

0.031 

0.526 

Next,  define  a current  velocity  u(z,t)  for  t > 0 and  R(t),  the  integrated  velocity 
vector,  such  that 


R(t)  = f°  u(z,t)  dz 

—oo 


(38) 
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Figure  9A.  Counterclockwise  spectrum  of  the  February  1974  current  record.  D,  f, 
and  SD  represent  diurnal  tidal  frequency , inertial  frequency,  and  semidiurnal  tidal 
frequency,  respectively. 


The  vector  R(t)  tends  to  decay  in  amplitude  under  the  influence  of  viscosity  and  to 
rotate  at  the  natural  frequency  of  the  system,  i.e. , the  inertial  frequency.  If  at 
a given  depth,  z,  uo(z)  leads  Rq,  then  for  t > 0 the  period  of  rotation  of  u(z,t) 
is  less  than  the  inertial  period  as  u(z,t)  tends  to  R(t);  if  Uq(z)  lags  behind  Rq, 
the  period  of  rotation  of  u(z,t)  is  greater  than  the  inertial  period  (Gonella, 

1971).  In  either  case,  the  inertial  period  is  approached  in  the  limit  at  each 
depth.  Gonella's  solution  (Gonella,  1971)  for  a wind  stress  imposed  at  to  = 0,  as 
well  as  for  a duration-limited  wind  stress,  predicts  the  same  approach  of  the  rota- 
tion period  to  the  inertial  period.  Figure  10,  adapted  from  Gonella  (1971),  illus- 
trates this  point.  Gonella  (1971)  noted  that  the  action  of  viscosity  explains  in 
part  peaks  in  his  spectrum  for  f - f i 2 to  4%,  where  f'  is  observed  frequency. 

The  same  mechanism  seems  to  be  responsible  for  the  spectral  peak  spanning  the  fre- 
quencies 0.041  to  0.043  CPH,  which  are  within  f = f + 6%. 

Modification  of  the  period  of  inertial  rotation  may  result  from  a number  of 
conditions,  among  which  are  mass  stratification  or  the  presence  of  horizontal  mean 
flow  (Gonella,  1971;  Perkins,  1972).  The  latter  may  create  a Doppler  shift  in  the 
observed  frequencies.  Perkins  (1972)  noted  a frequency  of  rotation  3 percent  above 
the  theoretical  f because  of  the  presence  of  a mean  current.  Gonella  (1971) 
observed  similar  frequency  shifting  in  three  examples  of  Inertial  motion.  Pollard's 
(1970)  wind-driven  model  for  a two-layered  ocean  predicts  frequencies  in  the  iner- 
tial current  in  the  range  f < f'  < 1.07  f.  The  fact  that  in  this  particular  study 
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Figure  9B.  Clockwise  spectrum  of  the  February  1974  current  record.  Inertial- 
diurnal  peak  is  indicated  at  top.  D,  f,  and  SD  represent  diurnal  tidal  frequency, 
inertial  frequency,  and  semidiurnal  tidal  frequency,  respectively. 


the  major  spectral  peak  occurs  at  a frequency  6 percent  above  inertial  is  not  sur- 
prising. 

Stratification  probably  contributes  to  some  extent  to  the  frequency  shifting. 
Stratification  of  waters  in  the  study  site  is  not  uncommon  because  of  land  runoff 
and  Mississippi  River  discharge.  This  stratification,  as  well  as  the  presence  of 
mean  flow,  operates  in  a rather  complex  fashion  that  is  difficult  to  define  quanti- 
tatively. For  this  reason,  shifting  of  the  inertial  frequency  is  usually  treated 
in  the  literature  only  in  a very  general  and  qualitative  manner  (see,  for  example, 
Perkins,  1972;  Gonella,  1971). 

Although  the  major  peak  at  0.043  CPH  may  be  explained  by  the  effects  of  mass 
stratification  and  Doppler  shifting  of  the  inertial  frequency,  the  reason  for  the 
four  minor  peaks  at  0.031,  0.035,  0.049,  and  0.053  CPH  is  not  as  clear.  On  the 
basis  of  available  data,  it  is  difficult  to  make  definitive  conclusions  as  to  their 
cause.  However,  it  does  seem  likely  that  they  may  be  a result  of  some  property  of 
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Figure  9C.  Total  spectrum  of  the  February  1974  current  record.  D,  f,  and  SD  repre- 
sent diurnal  tidal  frequency,  inertial  frequency,  and  semidiurnal  tidal  frequency, 
respectively. 


Figure  10.  Projection  of  Ekman  spiral  on  a horizontal  plane.  Rotation  period  for 
u(z,t)  above  Ro  is  less  than  inertial,  and  below  Rq  it  is  greater  than  inertial. 
(Adapted  from  Conella,  1971.) 
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the  wind  field.  Each  of  the  four  frequencies  has  an  associated  Cr  > 0.93,  suggest- 
ing strongly  polarized  clockwise  motion.  The  four  peaks  may,  therefore,  represent 
driving  by  a clockwise  rotation  of  a strong  wind.  Each  of  the  four  frequencies 
suggests  a separate  wind  event  and  strong  energy  at  a particular  frequency. 


E.  Simulation  of  Wind-Induced  Currents  with 
the  Pollard-Millard  Model 


Pollard  and  Millard  (1970)  were  able  to  account  for  wind-stress- induced 
inertial  currents  in  the  surface  layers  of  the  open  ocean  with  a mathematically 
simple  model.  They  assumed  a homogeneous  surface  layer  of  uniform  depth  on  an  f- 
plane.  Rather  than  treating  the  wind  stress  as  a surface  boundary  condition,  they 
assumed  it  to  be  instantaneously  and  uniformly  distributed  through  the  mixed  layer 
of  depth  Z,  i.e.,  as  a body  force.  This  particular  modeling  of  the  wind  stress  is 
justified  for  two  reasons  (Pollard,  1970): 

1,  The  depth  of  the  well-mixed  surface  layer  created  by  applying  a surface 
stress  increases  as  tl/3  (Rato  and  Phillips,  1969),  where  t is  time.  On  the  other 
hand,  the  time  scale  for  transmitting  momentum  through  the  well-mixed  layer  several 
tens  of  meters  deep  is  on  the  order  of  several  hours.  Therefore,  the  time  scale 
for  spreading  the  momentum  through  the  mixed  layer  is  short  compared  with  the  iner- 
tial period,  but  the  time  scale  of  the  deepening  of  the  layer  is  considerably 
longer  than  the  inertial  period.  Consider,  for  example,  a constant  wind  stress  of 
T = 4 dynes/cm^  applied  at  the  surface  of  a strongly  stratified  layer  of  Brunt- 
Vaisaia  frequency  N = lO”^  sec“l.  Kato  and  Phillips'  formula  states  that  the  mixed 
depth  Z(t)  is  given  by  Z(t)  = u*  (15t/n2)l/3^  where  u*  = (x/Py,)l/2  is  the  friction 
velocity,  t is  time,  and  p^j  is  water  density.  A wind  stress  of  4 dynes/cm^  would 
erode  such  a stratified  fluid  down  to  20  ra  in  less  than  2 hours.  On  the  other  hand, 
the  rate  at  which  Z(t)  increases  is  given  by  the  derivative  of  the  above  equation, 
dZ(t)/dt  = 0.82  u*/(tN)2/3.  After  6 hours,  the  constant  wind  stress  deepens  the 
layer  at  a rate  of  1.6  m/hr;  after  12  hours,  that  rate  decreases  to  1 m/hr.  In 
other  words,  the  time  required  for  the  bottom  of  the  surface  layer  to  "feel"  momen- 
tum applied  at  the  surface  is  at  least  an  order  of  magnitude  shorter  than  the  time 
for  the  depth  of  the  layer  bottom  to  increase  substantially.  Instantaneous  trans- 
mission of  surface  stress  through  the  surface  layer  is  therefore  a valid  assumption. 

2.  A weakly  stratified  layer  cannot  sustain  a strong  vertical  velocity  shear 
and  therefore  will  have  a tendency  to  have  momentum  uniformly  distributed  through 
it.  Homogeneous  distribution  of  momentum  from  the  surface,  where  stress  is  applied, 
to  a depth  Z is  thus  a valid  assumption. 

The  vertically  averaged,  linearized  equations  of  motion  in  the  mixed  layer 

are: 


F - ku 

u,  v = (u',v')  dz  (39) 

G - kv 

where  u'  and  v'  are  point  measurements  of  the  current  and  u and  v represent  their 
vertical  averages  and  F and  G are  the  forcing  functions  (in  this  case,  wind  stress) 
in  the  x and  y directions,  respectively.  The  F and  G time  series  are  here  defined 
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where  sin  0^  and  U^  cos  0^  are  the  1-hour  vector  averages  of  wind  velocity  com- 
ponents for  hours  i = 1,  2,  3,...n,  and  p„  are  the  density  of  air  and  water, 
respectiveiy , Z is  the  depth  of  the  mixed  water  layer,  and  Cy  is  the  drag  coeffi- 
cient (=  1.2  X 10”^  (Hsu,  197A)).  The  linear  damping  terms  -ku  and  -kv  introduce 
an  exponential  decay  factor  into  the  solution  of  equations  39  (see  "Introduction"). 
I'he  constant  k~^,  referred  to  as  the  e-folding  time,  describes  the  Lime  required  for 
the  rotary  current  to  be  dissipated  by  viscosity  to  e”^  of  its  original  value.  If 
k"^  = T,  where  T is  an  arbitrary  period  of  time,  then  after  a time  t = T w""  have 
the  current  speed  V(t)  related  to  the  initial  speed  by  V(t)  = Voe“(l/T)'I  = V^e'^ . 
Equations  39  were  solved  iteratively  by  Pollard  and  Millard  using  a time  step  of 
1 hour.  Considering  the  simplicity  of  the  model,  the  predicted  currents  they 
obtained  show  striking  similarity  to  those  observed. 

The  Pollard-Millard  model  was  used  to  generate  currents  from  Weather  Bureau 
wind  data  at  Boothville  for  28  January  to  4 February  and  winds  from  the  Chevron  rig 
for  4 through  28  February.  Profiles  of  temperature  and  salinity  near  the  moored 
meter  at  the  beginning  and  tlie  end  of  this  period  indicate  a fairly  well  mixed 
surface  layer  down  to  about  10  meters  for  the  current  record.  Several  values  of 
the  e-folding  time  were  tried  in  order  to  obtain  the  best  possible  correspondence 
between  observed  and  modeled  current  amplitudes. 

Keeping  in  mind  the  points  made  in  Section  1 1-C  about  properties  of  the  wind 
field  that  Initiate  or  destroy  inertial  motions,  a comparison  of  wind  (Fig.  ilA) 
with  the  observed  and  modeled  current  (Fig.  IIB)  reveals  several  examples  of  wind- 
generated oscillations.  The  best  example  of  this  phenomenon  is  the  very  regular 
sinusoidal  oscillation  of  the  u and  v components  of  the  observed  current  beginning 
on  9 February  and  continuing  for  nearly  a week.  The  model  is  highly  successful  in 
reproducing  these  observed  oscillations.  A significant  feature  of  the  modeled 
current  is  the  fact  that  it  is  in  phase  with  the  measured  current  immediately  fol- 
lowing a wind  event  but  gradually  drifts  out  of  phase.  Pollard  and  Millard  (1970) 
noted  that  the  dominant  frequencies  of  the  model  and  the  ocean  are  slightly  differ- 
ent, and  therefore  the  two  will  drift  out  of  phase  when  there  is  little  forcing. 
Furthermore,  the  damping  factor,  k,  may  change  significantly  with  time,  thereby 
altering  the  character  of  the  observed  oscillations.  As  a result,  the  model  can 
be  expected  to  give  the  best  correlation  with  the  observed  current  within  a few 
days  of  a given  wind  event. 

Figure  IIB  shows  the  comparison  between  the  simulated  current  with  2-day 
damping  (k  = 5.8  x 10“^  sec“^)  and  the  observed  current.  The  model  reproduces  with 
considerable  accuracy  nearly  all  the  major  features  of  the  measured  current  for  the 
period  4 through  28  February.  It  is  not  as  successful  over  the  period  28  January 
to  4 February  for  reasons  to  be  discussed  later. 

On  9 February  a moderately  strong  wind  rotating  nearly  inertially  in  phase 
with  the  current  continuously  added  momentum  to  the  current  for  a 24-hour  period, 
fhe  result  was  the  set  of  large-amplitude  oscillations  that  gradually  diminished  in 
strength  with  time.  Wind  cn-'ditions  following  the  wind  shift  on  the  9th  and  con- 
tinuing through  the  15th  were  very  conducive  to  maintaining  the  oscillatory  motion; 
that  is,  after  the  9th  and  up  to  the  15th  the  wind  was  weak  and  therefore  there  was 
little  additional  forcing  to  either  enhance  or  destroy  the  current  motion. 
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Figure  11.  A.  Wind  record  at  Chevron  platform.  B.  E-W  and  N-S  components  of  observed  (solid  line)  and 
simulated  (dotted  line)  currents.  F in  A indicates  approxinuite  time  of  frontal  passage.  Simulated  current 
from  28  January  to  the  afternoon  of  A February  was  generated  using  Weather  Bureau  wind  record  at  Boothville. 
From  4 to  28  February  wind  at  the  Chevron  platform  was  used.  Dates  are  for  January  and  February  1974. 


Without  knowledge  of  the  wind  characteristics  during  this  period,  the 
observed  oscillations  might  be  assumed  to  be  diurnal  tidal  motion.  In  fact,  fil- 
tering with  the  Doodson-Warburg  and  binomial  filters  seems  to  indicate  that  this  is 
indeed  the  case  (see  Fig.  5).  However,  the  possibility  that  these  are  purely 
tidally  driven  motions  must  be  ruled  out  because  tlie  tidal  record  at  Bayou  Rigaud 
(Fig.  6B) , as  well  as  the  NOS-predicted  tide  for  Caminada  Pass  (National  Ocean  Sur- 
vey, 1974)  (Fig.  6A) , indicates  relatively  small  and  partially  semidiurnal  tidal 
amplitude  for  tliis  period.  The  period  corresponds  to  the  time  of  equatorial  tides, 
when  tidal  amplitude  was  at  a minimum.  It  must  be  remembered  that  the  Pollard- 
Millard  model  is  a purely  wind-driven  model  and  therefore  it  cannot  produce  rotary 
tidal  motion.  The  degree  to  which  it  accurately  reproduces  the  observed  oscilla- 
tions of  9 through  15  February  is  strong  evidence  that  the  observed  rotary  motion 
was  wind  Induced. 

On  15  February  a frontal  passage  late  in  the  day  was  ."ccompanied  by  a 360° 
turn  in  the  wind  vector  over  approximately  a 6-hour  period.  Its  effect  was  the 
destruction  of  whatever  remained  of  the  rotary  motion  induced  by  the  9 February 
wind  event.  The  observed  current  on  the  16th  and  17th  showed  some  semidiurnal 
variation  that  was  not  apparent  in  the  simulated  current.  These  may  represent  an 
effect  of  either  the  semidiurnal  tide  or  turbulence. 

Anticyclonic  motion  was  created  by  a wind  shift  on  17  February  and  was  some- 
what enhanced  by  the  front  late  on  the  18th.  Although  this  front  initially  produced 
moderate  to  strong  currents,  the  oscillatory  amplitude  was  diminished  from  noon  on 
the  19th  through  midnight  on  the  20th.  The  reason  is  that  the  wind  was  very  strong 
and  nearly  unidirectional  for  about  one  inertial  period  on  the  19th.  Wind  during 
this  period  seems  to  have  removed  momentum  from  the  current  system. 

On  the  afternoon  of  21  February  a front  initially  produced  strong  amplitudes 
in  the  current  components  that  were  followed  by  undulations  of  near-inertial  period 
that  lasted  for  about  two  and  a half  inertial  periods.  These  undulations,  however, 
were  not  as  strong  in  amplitude  as  might  have  been  expected  if  we  examine  the 
strength  of  the  wind  stress.  As  before,  the  reason  appears  to  be  the  direction  of 
the  wind  on  the  21st  and  22nd.  The  wind  was  unidirectional  for  nearly  one  inertial 
period  beginning  on  the  21st  and  ending  on  the  22nd.  It  was  therefore  instrumental 
in  diminishing  the  oscillatory  amplitude. 

A word  of  caution  should  be  made  at  this  point  regarding  the  spacing  of  wind 
events  and  the  oscillations  they  produce  in  the  current  components.  As  one 
approaches  lower  latitudes,  the  inertial  period  becomes  Increasingly  long.  At  the 
latitude  of  the  study  site,  this  period  is  slightly  longer  than  a day.  If  wind 
events  occur  one  or  two  inertial  periods  apart,  as  is  the  case  in  the  instance  just 
mentioned,  it  becomes  increasingly  difficult  to  observe  inertial  motions.  In  the 
strict  sense,  inertial  currents,  once  induced,  are  acted  upon  only  by  the  deflecting 
force,  Coriolis  force.  In  real-world  situations,  inertial  rotary  currents,  once 
Induced,  are  acted  upon  by  both  Internal  and  external  forces,  which  may  make  their 
observation  difficult.  For  example,  if  a front  creates  what  appears  to  be  inertial 
rotary  motion  and  is  followed  within  one  or  two  inertial  periods  by  another  front 
and  accompanying  clockwise  wind  rotation,  the  current  motion  observed  between 
frontal  passages  may  not  be  Inertial  rotation  but  rather  currents  that  are  driven 
by  Che  ancicyclonically  rotating  wind  vector  associated  with  the  second  approaching 
front.  As  a result,  characterization  of  the  oscillations  between  times  of  the 
fronts  of  21  and  24  February,  Cor  example,  may  be  somewhat  problematical  because  of 
the  clockwise  rotation  of  the  wind  on  the  21st,  which  was  followed  by  another  rota- 
tion on  the  22nd  and  a third  on  the  24th. 

Strong-amplitude  oscillations  appearing  in  the  observed  current  toward  the  end 
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of  Che  record  were  induced  by  Che  fronc  of  24  February.  Winds  during  chis  period 
were  generally  weak  Co  moderace  and  shifced  very  slowly  clockwise.  The  degree  Co 
which  Che  model  is  able  Co  simulaCe  Che  observacions  is  only  fair.  Fspecially  sur- 
prising is  Che  magnicude  of  Che  observed  currenc  on  26  February,  ciie  sCrongesc  cur- 
renc  of  Che  encire  scudy  period.  Alchough  die  model  does  produce  significanC  undu- 
lacions,  ic  falls  far  shore  of  creating  che  peak  of  26  February.  As  was  die  case 
earlier  in  che  record,  cidally  driven  currencs  must  be  ruled  ouC  as  a primary  mecha- 
nism in  producing  the  peak.  The  tide  record  at  Bayou  Kigaud  shows  relatively  small 
cidal  range  for  25  to  27  February,  as  is  the  case  with  the  NOS-predicced  tide  levels. 
This  particular  period  was  the  time  of  equatorial  tides  and  therefore  tidal  effects 
were  minimal.  The  failure  of  the  model  to  reproduce  die  strong  magnitude  of  the 
rotary  current  of  the  26th  may  possibly  be  explained  by  a significantly  reduced 
friction  coefficient,  k,  for  about  a day,  but  there  are  no  data  available  to  sup- 
port this  conjecture. 

One  final  event  worth  mentioning  is  the  front  of  3 February.  Unfortunately, 
recording  of  wind  data  at  the  Chevron  rig  was  not  begun  until  the  afternoon  of  the 
4ch.  The  observed  current  components  following  the  front  show  fairly  well  behaved 
strong  oscillations,  and  it  was  therefore  desirable  to  model  them.  Using  1-hour 
values  of  wind  velocity  collected  by  the  Weather  Bureau  at  Boothville,  simulated 
currents  were  generated  for  the  period  before  4 February  not  covered  by  the  Chevron 
wind  data.  A linear  regression  of  che  available  Chevron  wind  data  was  performed  on 
the  boothville  wind  data  to  compensate  for  differences  in  wind  magnitude  over  land 
and  over  water.  The  results  at  1500  hours  on  4 February  of  modeling  with  the 
regressed  wind  data  were  used  as  initial  current  amplitudes  for  the  model  using  the 
Chev  on  wind  record.  This  procedure  produced  only  marginally  accurate  results. 
Although  undulations  do  occur  in  the  simulated  components  beginning  on  3 February, 
they  do  not  adequately  describe  the  well-defined  oscillations  of  the  observed  com- 
ponents. In  addition,  currents  before  4 February  which  are  generated  with  the 
Boothville  wind  as  input  are  a rather  poor  simulation.  The  problem  stems  from  a 
number  of  causes,  among  which  are  (1)  winds  over  land  are  deflected  more  from  their 
geostrophic  values  than  those  over  water  and  therefore  wind  direction  at  Boothville 
is  probably  not  a good  representation  of  direction  over  the  current  meter;  (2)  dis- 
tance between  Boothville  and  the  moored  current  meter  is  Coo  great  to  allow  an 
adequate  representation  of  wind  conditions  at  the  study  site.  Despite  these  diffi- 
culties the  model  indicates  that  the  undulations  in  the  current  components  observed 
for  4 or  5 days  after  the  front  of  3 February  were  wind- induced  inertial  oscillations. 


III.  SUMMVRY  AND  CONCLUSIONS 

Results  of  the  various  analyses  performed  on  the  current  velocity  time  series 
suggest  strongly  that  the  rotary  motion  observed  is  locally  wind-induced  inertial 
oscillation  and  is  not  the  result  of  astronomical  tides.  Evidence  to  support  this 
conclusions  is  as  follows: 

1.  The  observed  currents,  plotted  as  a progressive  vector  diagram,  clearly 
show  clockwise  rotation  following  significant  wind  shifts. 

2.  Conventional  analysis  using  filters  to  extract  tidal  currents  from  the 
observed  time  series  produces  current  hodographs  that  are  totally  unrelated  to  tlie 
predicted  tide  or  to  the  measured  tidal  elevations.  The  size  of  the  hodographs  is, 
however,  closely  related  to  wind  shifting  usually  associated  with  frontal  passage. 

3.  Complex  demodulation  of  the  current  time  series  indicates  a large- 
magnitude  signal  near  the  inertial  frequency  that  is  closely  associated  with  passing 
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fronts.  The  major  peaks  have  a magnitude  of  30-40  cm/sec,  whereas  the  background 
magnitude  is  0-10  cm/sec  (see  Fig.  7).  For  this  particular  study  period,  the  latter 
is  probably  a good  measure  of  the  magnitude  of  rotary  currents  not  induced  by  wind, 
that  is,  tidal  currents.  Furthermore,  the  phase  of  the  demodulated  current  does  not 
show  a slow  uniform  decrease,  as  would  be  expected  for  tidal  currents.  Instead,  it 
demonstrates  the  transient  nature  of  locally  wind-induced  inertial  oscillations. 

The  phase  achieves  a constant  value  with  each  magnitude  peak  and  is  different  for 
each  of  those  peaks. 

4.  Spectral  analysis  indicates  a strong  peak  in  the  inertial-diurnal  fre- 
quency band.  The  record  is  not  long  enough  to  resolve  the  inertial  and  diurnal 
frequencies,  but  the  peak  appears  to  be  too  strong  to  be  generated  by  tides  alone. 

5.  The  Pollard-Millard  model,  which  uses  wind  stress  as  the  sole  input,  gives 
an  excellent  reproduction  of  the  observed  current  components.  Oscillations  are  pro- 
duced in  the  modeled  current  by  shifting  winds.  Particularly  significant  is  the 
agreement  between  modeled  and  observed  current  components  following  the  wind  shift 
of  9 February.  The  generally  good  agreement  in  magnitude  of  the  modeled  and 
observed  currents  indicates  that  wind  forcing  is  considerably  more  important  than 
tidal  forcing  in  producing  the  observed  rotary  motion. 

The  generation  of  inertial  oscillations  by  wind  is  dependent  on  two  properties 
of  the  wind  field:  its  duration  and  its  rotation  rate.  A unidirectional  wind 
generates  inertial  rotary  motion  during  the  first  half  of  an  inertial  period  and 
diminishes  it  during  the  second  half.  A wind  vector  rotating  at  the  inertial  fre- 
quency generates  an  inertial  current  and  continuously  increases  its  magnitude.  The 
magnitude  of  inertial  oscillations  observed  following  a particular  wind  event  is 
also  dependent  on  clockwise  rotary  motion  near  the  inertial  frequency  existing 
prior  to  the  wind  event.  A particular  observation  of  "inertial"  motion  may  there- 
fore represent  an  enhancement  of  inertial  oscillations  or  diurnal  tidal  currents 
existing  prior  to  a given  wind  event.  Not  only  is  the  wind  instrumental  in  gener- 
ating inertial  oscillations,  it  can  also  diminish  or  destroy  them.  If  the  wind 
vector  has  a component  opposite  the  current  vector,  it  will  remove  momentum  from 
the  current  field.  If  the  situation  persists,  clockwise  rotary  motion  of  the  cur- 
rent is  totally  eliminated. 

Inertial  motion  as  depicted  by  the  demodulated  current  (Figs.  7 and  8)  is  con- 
taminated by  diurnal  tidal  motion.  Inasmuch  as  the  inertial  and  diurnal  tidal  fre- 
quencies for  the  Louisiana  coast  are  so  close,  complex  demodulation  is  unsuccessful 
in  separating  the  two.  However,  demodulation  clearly  shows  that  wind  has  a dominant 
effect  in  generating  rotary  currents  in  the  inertial-diurnal  frequency  band.  This 
effect  is  apparent  when  we  examine  the  magnitude  and  phase  of  the  demodulated  cur- 
rent in  the  presence  and  absence  of  wind  events.  The  conclusion  may  be  further 
substantiated  by  examining  the  tide  record  for  the  region  during  the  study  period. 

The  closeness  of  the  inertial  and  diurnal  tidal  frequencies  also  presents 
problems  in  the  spectral  analysis.  The  32-day  record  examined  here  is  not  long 
enough  to  resolve  the  two  frequencies  as  distinct  spectrum  lines.  However,  even 
with  considerably  better  resolution,  formidable  problems  would  be  encountered.  As 
pointed  out  in  Section  II-D,  inertial  motion  is  observed  not  at  a single  frequency 
but  over  a band  of  frequencies  near  the  theoretical  inertial  frequency.  This 
"smearing"  of  current  energy  over  a frequency  band,  as  pointed  out  by  Gonella  (1971), 
results  from  the  observed  rotary  current  at  a given  depth  approaching  the  inertial 
frequency.  Therefore,  energy  associated  with  a diurnal  tidal  current  would  not 
appear  as  a distinct  peak  but  would  be  contaminated  by  energy  associated  with  the 
inertial  current.  Further  difficulty  would  be  introduced  by  shifting  of  frequencies 
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caused  by  the  presence  of  mean  flow  and  mass  stratification. 

Wind-induced  inertial  oscillations  of  the  surface  layer  along  the  Louisiana 
coast  appear  to  be  a fairly  common  phenomenon  if  the  record  examined  here  is  any 
indication.  Because  atmospheric  fronts  are  most  common  and  strongest  along  the 
Louisiana  coast  during  the  winter,  early  spring,  and  late  fall  (Leipper,  195A), 
these  periods  can  be  expected  to  produce  inertial  rotary  motion  frequently.  The 
frequency  of  frontal  passage  through  tlie  region  during  this  p riod  is  one  every 
3 to  6 days.  Thus,  inertial  currents  are  probably  very  commonplace  for  part  of  the 
year.  For  example,  the  32-day  winter  record  examined  here  contains  inertial  cur- 
rents most  of  the  time. 

On  the  other  hand,  atmospheric  fronts  decrease  in  frequency  and  strength 
through  the  late  spring,  sununer,  and  early  fall,  when  the  Bermuda  high-pressure 
system  dominates  the  weather  along  coastal  Louisiana.  Therefore,  inertial  currents 
are  much  less  likely  to  be  observed  during  this  period.  As  an  example,  the  39-day 
record  beginning  25  April  and  ending  2 June  contains  only  one  distinct  case  of 
inertial  motion. 
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Two  current  vector  time  series  obtained  in  the  Mississippi  Bight  exhibit  clock- 
wise polarized  currents  of  near-inertial  frequency  that  are  closely  associated  with 
shifting  winds.  Because  of  the  closeness  of  the  local  inertial  period  and  the  diurnal 
tidal  period,  it  is  difficult  at  first  glance  to  determine  the  true  nature  of  the 
observed  rotary  currents.  However,  complex  demodulation  at  the  inertial  frequency 
reveals  a strong  signal  accompanying  wind  shifts  that  are  usually  associated  with  the 
passage  of  atmospheric  fronts.  , 

Spectral  analysis  for  clockwise  and  counterclockwise  frequencies  indicates  a 
highly  energetic  peak  in  the  inertial-diurnal  frequency  band  for  the  clockwise  spectrum 
The  rotary  coefficient  computed  from  the  autospectra  and  quadrature  spectrum  of  the 
jvector  components  gives  Cr  > 0.9  in  the  vicinity  of  the  inertial-diurnal  frequency  band. 

A model  using  wind  stress  as  a forcing  function  is  highly  effective  in  repro- 
Iducing  sinusoidal  oscillations  seen  in  the  observed  current.  These  oscillations  occur 
|in  conjunction  with  shifts  in  the  wind  direction.  (11) 

Because  of  the  close  association  of  the  near-inertial  oscillations  with  local 

tlnd  effects,  it  is  concluded  that  inertial  currents  are  locally  induced  by  wind  stress 
urthermore,  wind  stress  not  only  Initiates  the  rotary  currents  but  is  also  highly 
ffectlve  in  destroying  tliem.  (U) 
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